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The nonlinear dynamics of wave envelopes modulated in 2+1 dimensions is 
considered for two systems in plasma physics: (i) Langmuir pulses and (ii) 
intense (but weakly relativistic) electromagnetic (EM) pulses. Using singular 
perturbation techniques applied to an envelope approximation, both problems 
are reduced to the two-dimensional nonlinear Schrödinger (2DNLS) system, 
which describes the dynamics of two coupled slowly varying potentials. The 
general 2DNLS system exhibits a rich variety of phenomena, including 
enhanced (compared with ‘longitudinal’ propagation) modulational stability 
and (1D) soliton formation; decay of 1D solitons over long time scales; self- 
focusing regimes (determined by a virial-type condition); as well as integrability 
and 2D solitons. Applying our recent results on the 2DNLS system, we 
determine which of these phenomena can actually occur here and compute the 
parameter regimes. (i) The 2DNLS system for the Zakharov equations is 
modulationally unstable for all parameter values. It also has an integrable 
sector and a self-focusing regime. (ii) The 2DNLS system describes coupled 
‘longitudinal’ and ‘transverse’ modulations of linearly or circularly polarized 
EM pulses propagating through a warm unmagnetized two-component neutral 
plasma with arbitrary masses (i.e. electron—positron or electron-ion). The pulse 
can accelerate particles to weakly (but not fully) relativistic velocities; 
relativistic, ponderomotive and harmonic effects all contribute to the nonlinear 
terms. The resulting 2DNLS system does not admit a self-focusing regime. 
Parameter values leading to an integrable case (the so-called 'Davey-- 
Stewartson I’ equations, which admit 2D soliton solutions) are computed; 
however, the required values would not be attainable in a laboratory or 
astrophysical setting. None the less, the existence of new nonlinear mod- 
ulational instabilities associated with the second spatial degree of freedom 
already represents an important potential limitation on any (1+ 1)-dimensional 
approach to nonlinear evolution and modulational instability of plasma EM 
waves. 
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1. Introduction and Model 


In a series of paper (Kates & Kaup 1989a, b, 1991, 1992, 1993; hereinafter 
referred to as papers I-V) we have seen that problems involving the intense 
propagation of plane-fronted EM pulses are more subtle than they appear at 
first glance, and require careful application of singular perturbation methods. 
The evolution of a modulated pulse is determined in general by the combined 
effects of linear dispersion, relativistic nonlinearity, nonlinear harmonic 
interactions and ponderomotive nonlinearity. (Harmonie interactions vanish 
identically for circular polarization.) As a consequence, all of these effects may 
be significant in computing criteria for nonlinear modulational instabilities 
(Zakharov & Rubenchik 1972) and possible soliton formation. 

Unfortunately, (1+1)-dimensional (‘1D’) treatments of modulated EM 
pulses provide only part of the information needed for applications: in addition 
to modulation along the propagation direction, transverse modulations may be 
excited, i.e. the wave fronts may be slightly curved. This area of investigation 
has taken on added importance recently in view of recent experiments 
involving ‘channeled propagation’ (Borisov et al. 1992a). 

In those systems for which plane-fronted pulses satisfy the scalar 1D 
nonlinear Schrédinger (NLS) equation 


iq, q,, t lg q = 0 (1) 
(see papers I, II and V) one might have supposed that generalization from 1D 
to 2D modulations of EM waves could have been accomplished simply by 
including transverse dispersion, i.e. by replacing the 1D dispersion term q,, in 
the NLS equation (1) with a 2D dispersion term: 


iq, + (Dzz 0; + 2D,y 0, Oy + Dy, %) a lal’ q = 0, (2a) 
where Dy = JURE (25) 
V 0k; Ok, 


A system of the form (2) is in fact obtained for a circularly polarized beam 
with modulations in both transverse directions but vanishing longitudinal 
modulations (Spatschek 1977). This case has been called the ‘thin-beam’ 
approximation. 

Generally speaking, however, a proper description of 2D modulations of EM 
pulses involves two independent cowpled potential functions. This point was 
recognized by Karpman (1990). In particular, as we shall shortly see, (2) will not 
provide the correct 2D evolution equation for a warm plasma if both 
longitudinal and transverse modulations are present, because fluid modes 
associated with longitudinal and transverse ponderomotive forces are coupled. 

With this motivation, we recently studied (Kates & Kaup, 1994; hereinafter 
referred to as paper 0) a general class of two-dimensional nonlinear Schrödinger 
(2DNLS) systems, which can all be expressed in the ‘canonical’ form 


10,9+0,9 = pa, (3a) 

0,p = Os(q*q), (36) 

where each of the operators O, in (3) is a second-order linear symmetric 
(dispersion) operator of the form 

0, = DO? 02 - 2D™ 0,0, +D 02 (3c) 


yy "y" 
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The 2DNLS system (3) subsumes many special cases of interest, as described in 
paper 0. These include the ‘DR’ system (Djordjevic & Redekopp 1977; 
Ablowitz & Segur 1979), first studied in the context of water waves, and of 
considerable interest in its own right: 


iA, +AA g+ LA, = XA*A? t y, A®,, (4a) 
a0,,--0,, — —B(A*A),. (4b) 


Any (1+ 1)-dimensional restriction of (3) (or (4)) satisfies (1). 

The 2DNLS system exhibits a rich variety of phenomena in various regimes, 
including enhanced modulational instability, self-focusing, integrability and 2D 
soliton formation. The regimes of self-focusing and modulational stability are 
mutually exclusive, as are the conditions for integrability and self-focusing. 

As shown in paper 0 (see (0.30)+ the 2DNLS system is modulationally stable 
(for given operators, O,, n = 1,2,3) provided that 


g mg oo 0 VK, (ba) 
where K is the wavenumber of the modulation, and where 


See (3) for the definition of the Df?? in terms of the O,. 

If (5a) is not satisfied for some K then the system is modulationally unstable 
to modulations in the direction K (see paper 0) for sufficiently small K (= |K]). 
The 1D projection (defined by assuming spatial dependence on y = K.x) 
satisfies the NLS equation (1) with the upper (modulationally unstable) sign, 
and thus admits 1D soliton solutions. 

The self-focusing regime of (3) is characterized by the eigenvalues of the 
Di: O, and O, must both be elliptic, and for that choice of variables for which 
Df? and D both have positive signature, DP must have at least one negative 
eigenvalue. 

The integrable cases of the DR system, which are known as the 
Davey-Stewartson equations (DS-I and DS-II), as well as all other integrable 
cases of (3), may be obtained from the following conditions: first, 


at least one of O, and O, must not be elliptic; (6a) 
secondly, D? = xD? (65) 
must be satisfied for some constant «; thirdly 
EDP) a) = 0; (6c) 
ij 
and fourthly det(D™) +0 for each of n = 1,2,3. (6d) 


The cases admitting 2D solitons were discussed in paper 0. 

Consider a system that reduces to the (1D) NLS equation (1) in some 1D 
subspace (e.g. modulations in the z direction). If this ID reduction is 
modulationally stable, the system may still be unstable with respect to 
modulations in a different direction. In particular, it may self-focus or, if it 
happens to be integrable, form 2D solitons. Similarly, if a 1D reduction is 
modulationally unstable and forms 1D solitons, these may decay owing to 
transverse perturbations. These properties of (3) pose a potential limitation on 
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the validity of conclusions drawn from a 1D analysis of any intrinsically multi- 
dimensional, NLS-type system. 

In $2, we show that the canonical form (3) of the 2DNLS system describes the 
Zakharov equations (Zakharov 1972) in an envelope approximation. Although 
the Zakharov equations are electrostatic in nature, they are still an excellent 
model to use for illustrating the methods used in $3. Furthermore, this will 
provide an example of a non-trivial 2DNLS system in canonical form (3) that 
does not always reduce to the DR form (4). In this approximation the Zakharov 
equations are always modulationally unstable. They contain an integrable case 
as well as a self-focusing regime. 

In paper I we considered 1D modulations of a plane-fronted pulse with 
arbitrary polarization, propagating in an unmagnetized (possible warm) 
plasma. In $3 we consider the same model for the plasma, but allow 2D 
modulations. As in papers I, II and V, e = (e/mcw)|#| is a dimensionless 
measure, both of the typical velocity of electrons excited by the wave and of the 
slowness of the modulations, but here the ‘modulations’ depend on both y and 
z, where Z is the propagation direction. We give the explicit computation for 
linear polarization along X (circular polarization may be treated similarly). 
Consequences are discussed in §4: In the present case of an unmagnetized 
plasma it turns out that the effects of finite temperature are crucial in 
determining the qualitative behaviour. In particular, for an electron—positron 
plasma the frequency range for modulational instability including both 
transverse and longitudinal modes is considerably broader than it was for 
purely longitudinal modulations (paper I). 

We emphasize that in a self-consistent approximation procedure the effects 
of linear dispersion, relativistic nonlinearity, nonlinear harmonic interactions 
and ponderomotive nonlinearity will generally enter at the same order with 
respect to e, just as in the ({+1)-dimensional case. In the treatment of 
Karpman & Washimi (1977) relativistic effects were not included (see Shukla & 
Stenflo 1984). In the treatment of Borisov et al. (19925) the approximation 
p, = A stated below their equation (8) takes into account relativistic corrections 
(compare our (5) below) but not pondermotive forces (see (32)-(42) below or 
(I. 20) in the (1+ 1)-dimensional case). If one is interested in electron—positron 
plasmas, ponderomotive effects are at least as large as the other nonlinear 
effects, and must always be included. Now, in 1D systems, nonlinear 
ponderomotive effects in an ion-electron plasma are usually down by a factor of 
m/M compared with relativistic and harmonic effects. However, in the 
coefficient f in (45), which gives the strength of the coupling to the potential 
field ©, (transverse) ponderomotive forces contribute to the leading terms. 
Whether or not this coupling has a short enough time scale to be qualitatively 
important will of course depend on the details of the particular physical 
situation. 


2. Reduction of the Zakharov equations to the 2DNLS system 


The Zakharov equations (Zakharov, 1972) arose in models of the multi- 
dimensional interaction of acoustic waves with an electrostatic (Langmuir) 
plasma wave. Although they involve electrostatic and not electromagnetic 
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waves, they serve as a simple example of the application of the previous results. 
In fact, we can show that there is an other integrable sector of the Zakharov 
equations, which has the so-called DS-I form. Relativistic generalizations of the 
Zakharov equations have been discussed by Berezhiani, Tsintsadze & 
Tskhakaya (1980). As we shall note shortly, even the relativistic generalizations 
may be analysed by the same techniques. 

After normalization, the classical Zakharov equations for one linearly 
polarized mode may be written as 


0; P —c; VP = V*(Q*Q), (7a) 
i0,Q4- V1Q = PQ, (7b) 


where c, is proportional to the sonic velocity, and damping is not considered. 
Equations (7) exhibit collapse and other phenomena associated with non- 
integrable equations. If there exists a limit for which the 6? P term may simply 
be neglected, (7) evidently reduces to a pure (multi-dimensional) nonlinear 
Schródinger equation (2). However, as shown in paper 0, for integrability an 
auxiliary potential must be present. This may be achieved by a different 
limiting procedure: 

One expands P and Q in terms of a small parameter e, taking the leading term 
Q, in Q to be of first order and the leading term P, in P to be of second order 
in e. One then introduces slow co-ordinates 


X=ex, Y=ey, Tzd, reet (8) 
and assumes an envelope approximation for Q, of the form 
Qi = get(kr—ot) (9) 


(For convenience, the wave vector is chosen without loss of generality to be in 
the x direction. There is still slow dependence in both spatial directions.) Both 
P, and a are assumed to depend only on the slow variables. The first-order terms 


yield a dispersion relation 
9 E, (10) 


while elimination of secular terms at second order leads to the group velocity 
condition 


(0p T v,0y)a = 0, (110) 
do 
where (=a, = 2k. (115) 


The third-order secular condition on (7b) yields 
ia, +dyytayy = Pa, (12a) 


while the leading (third-order) term in (72) is 


v? 0% P,—c2(0% + 02) P, = (0x + df) a*a. (125) 

Equations (12) are of the canonical 2DNLS form (3), with 
0, = 0 0, (13) 
0, = (v5 — C5) Oy — 0s Oy, (14) 


0, = 0% + 05 =0,, (15) 
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The integrability conditions (6) then reduce to the requirement 
v? = 4k? = 26, (16) 


which gives the DS-I form. Thus we see that a sector exists for which the 
Zakharov equations are integrable in 2--1 dimensions: namely where the 
Langmuir wave is a small-amplitude almost-plane wave with group velocity 
near 2%,. Incidentally, this computation illustrates the utility of expressing 
integrability conditions for the ZDNLS system in a co-ordinate-free canonical 
form: this integrable sector of the Zakharov equations appears to have been 
unknown until now. 

Similar results hold for the relativistic generalizations of the Zakharov 
equations. Using the form given by Berezhiani ef al. (1980), if (7b) is replaced 
by 


10,Q+V?Q = PQ—TQ*Q* (10^) 
then the above procedure would also lead to the 2DNLS system but with 
0, = 0% + 0} — l'O,, (15’) 


In this case, O, need no longer be positive-definite. 

We can also use the canonical 2DNLS form (3) to study modulational 
instability by using (13)-(15) in (5). Since both Z2? and Z9? are positive- 
definite, modulational stability requires 


9290-0 (17) 


which, in view of (14), is always violated for some K. This result appears to be 
consistent with numerical investigations (Goldman 1984). 

Equation (12) satisfy the self-focusing conditions of paper 0, which are of the 
‘virial’ type, if 

và = 4k? < e. (18) 

This result may also shed some light on observations, and would seem to be 
consistent with the existence of self-similar solutions of the full equations 
exhibiting collapse (Goldman 1984). 


3. Reduction of EM-pulse propagation in a warm plasma to the 
2DNLS system 


In light of the preceding results, we now return to the problem considered in 
paper I: the long-term evolution of a modulated pulse propagating in an 
unmagnetized (possibly warm) plasma. We follow the notation of that paper 
unless otherwise stated. We take the z axis in the direction of propagation 
(k = zk) as in paper I, but permit slow modulations along both the z and y 
directions. In the following computation we assume linear polarization in the x 
direction. 

The system (L4)-(I:6) to be solved (Maxwell’s equations and two sets of fluid 
equations) describes a fully ionized plasma composed of singly ionized atoms (of 
charge +e and mass M) and electrons (of charge —e and mass m); the plasma 
is modelled as a perfect fluid with constant but possibly distinct electron and 
ion temperatures. No approximation is made with respect to m/M unless 
explicitly stated. (Upon setting M = m, we recover the electron—positron case.) 
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The background plasma is taken to be unmagnetized and neutral. (There will of 
course be deviations from neutrality due to nonlinear effects.) Units are chosen 
such that c = 1 unless otherwise stated. 

Mathematically, the intended generalization of paper I may be stated as 
follows: we study the evolution of appropriate initial data given in the plane 
—00 <z « 0, —% « y « oo for the functions h, v, N, V, A and ¢ (the electron 
density and velocity, the ion density and velocity, the vector potential and the 
electric potential respectively). We are interested in solutions that satisfy the 
original equations up to appropriate order (as explained in papers I-V) and that 
correspond to slowly modulated almost-plane waves. We assume that well- 
posed initial data can be given that evolve in this way. All results then follow 
from substitution of the expansions into the field equations and collection of 
terms order-by-order. Some of the equations are to be solved directly, while 
others yield integrability conditions. 

Of course, there will be regions of parameter space where three-wave resonant 
interactions could occur. These regions will usually be in a small frequency 
range near the plasma frequency. When they can occur, naturally they will tend 
to dominate over the four-wave interactions that give rise to the family of 
nonlinear Schródinger equations. In particular, Rizzato & Chian (1992) have 
recently studied the interaction of two oppositely propagating Langmuir waves 
with a transverse electromagnetic wave via the oscillating two-stream 
instability. These interactions can be very important in the regions of 
parameter space where they occur. However, here we shall assume that we are 
away from such regions of parameter space and shall restrict our attention to 
the four — wave interactions of the nonlinear Schródinger family of equations. 

Asymptotic expansions are assumed according to (L.7) and (I.8), where 
^, = const, Na =n, and all other zeroth-order quantities vanish. The slow 
variables are 

Zze, Y=, T=e, r=et. (19) 
With this scaling, the effects of linear dispersion (with respect to both z and y), 
relativistic corrections, harmonic interactions and ponderomotive forces will 
enter at the same (third) order of our approximation scheme. 

The linearized theory is given by (I.11)-(I.13), with 


A, = ae) 4 c.c, (20) 
and a in the (x,y) plane. Here we assume 
a — aX. (21) 


As in papers I-V, a will now be permitted to depend on the above slow co- 
ordinates. We collect like powers of e in all equations and enforce the non- 
secularity conditions of paper I (see (1.15) and (I.16)). As usual, the first-order 
terms are identical to the linear theory, while slow derivatives in a will enter the 
higher-order equations. Additional slowly varying quantities will be introduced 
as required. In particular, terms will be classified as ‘DC’ (0), ‘fundamental’ (1), 
‘harmonic’ (2), etc., as in papers I-V. For example, the second-order electron 
density n, is decomposed as 


Ny = NP + (nD) éD + o o.) + (ni? eik- + e o), (22) 


with analogous notation for other variables. Note that coefficients such as n? 
generally depend on the slow variables. 
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As before, we may take 
n® = NO = 6 — 0, 
a ey (23) 
v® = v(9 = AM — 0. 


without loss of generality. 
Elimination of secular terms at second order leads to the group velocity 
condition 


(05-9, 07) a = 0, (24a) 
do k 
where, as before, d (245) 


The absence of terms involving dy is simply due to the definition of z as the 
direction of propagation. 
The harmonic terms at second order may be shown to be given by 
2,2 
f nPHRND == Ra, (252) 


(2) = 2 
ny”  4mf,—mMo, 


NP ^ AMF,—mMo* e 
f, = m(w* — e$ ES, (25c) 
F, = M(w*—CkE?). (25d) 


These terms agree with those of paper I (given for vanishing pressure (I.21)). 
This is to be expected, since they depend only algebraically on a. 

The first-order DC parts of the fluid equations imply, with the second-order 
DC part of Poisson's equation (I.65), 


np = NY, (26) 
and with the second-order DC part of Maxwell’s equation (I.6c), 

V9 = v(9 (27) 

AQ = 0. (28) 


However, determination of the second-order DC quantities requires evaluation 
of the third-order fluid equations, and will be discussed shortly. 
One can again take 


pP =NY =n =0, AM =0 (29) 


without loss of generality. The relativistic corrections again produce third-order 
perturbation of the velocities of the form 


NE 3 

v® = ms a?a*X, (30a) 
3 

VO = zou. (305) 


(Note again that these derive from terms depending algebraically on a.) 
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From the fundamental part of the z component of the V x B equation at third 
order, we obtain the secular condition 


(0) 
Qiwd, at (102, +0 a} —d})a—wh a 
0 


a L NP 1 Submit eM sg yay 
050 e Mn, m+M- 28 mi4M1* a 0. (31) 


We now come to the second-order DC terms (n{®, N® etc.): in contrast to 
paper I, these may not all be expressed algebraically in terms of a (except in the 
case of vanishing pressure). We first introduce the abbreviations 


q = n$? [ny, (32) 
V, =ð, 4 26,. (33) 


The third-order DC parts of (I.45) and (I.55) may be added to eliminate the 
second-order DC potential, yielding 


(m 4- M) 0, v9 +e? (5; Vi(a*a) + (mc? - MC?) V, q = 0. (34) 

From the continuity equation (I.4a) at the same order, we have 
O5 qt V,. v — 0. (35) 
We now let v? = (0; y) f+ wz. (36) 


Note that v now has both y and z components, in contrast to (1.206). Like a, 
the quantities q, V, and w are functions of the slow variables. 

Since we are interested in solutions driven by a, we can eliminate T 
derivatives using 


67a = — v; za, (37a) 
Org = —v,07q, (37b) 
pw = —v, 02%, (37c) 
ôr y = —v,0z y, (37d) 
where (37a) is just (24a) and v, is given in (245). It is convenient to define 
C= mee (38) 


The z component of (34) may then be integrated to obtain 


2 
: - aa* +o, (39) 


g Ug 


while integration of the continuity equation (35) yields 


q= wt oy (40) 


Ug 


In each of the above integrations we have assumed that the various quantities 
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vanish at the front of the pulse and that the pulse is entering a quiescent region. 
Substitution of (39) into (40) then gives 


(net ay p). (41) 


Inserting the Z derivative of (37d) and the Y derivative of (41) into the y 
component of (34) then yields 


1 


=La 
v — o? 


ev 
(vg — v3) Vzz—VviVvy = zy mv (42) 


This coupled to the secular condition (31), which may be rewritten as 


; (95 C 2 "CO 
2iwa, +l Azz te Gy y t (Cg t 3Cg t Cp)a*a Bayly =0, (43) 
g 8 


where, as derived in paper I, Cy, Cp and C, are given by 


nnm eu? mM 1 F 1 
a 2€ M+m\m3(w?— kc?) | M3 (o? — kC?) 


o? mM 1 1 
urs kj Mlo- | | (14) 


. 2 0$ MM mM X aes He ER 
with D,z1— 4 M m|m(o? -di MoO) 
ew? m^? +M’ 
and Cr= 2c mM? on 
2 
Cp = p EEP E em) 


è Mm M(vi—Ci)m(v;—cj) 


The system of equations (42) and (43) may be recognized as the DR form of 
the 2DNLS system (0.3) if we identify 6 = y, A =a, = Y and y = Z and 
redefine 7 of (0.3) as 7/2w. The constant appearing in (4) are then given by 


x 7 —(Cy T 30g Cy), (47a) 
e zn 
x ~ mM(v? — v2) vy’ ei) 
A=1, (410) 
"má 
Ho = P (41d) 
Xi = vi—v : 
v? 
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If the computations of this section are repeated for the case of circular 
polarization, the system may again be reduced to the form (4), but (47a, b) are 
replaced by 


XxX7-—2(Cp-2Cg) (circular polarization), (47 a) 
f=- Ze (circular polarization) (47b) 
— mM (v2 —v?) P 


In terms of the canonical form (3), we have for both linear and circular 
polarization 


Lp s ô — vs J 
i rrt ( i-a a) (48a) 
Xv; —v$) zi Xs a2 
0, = "T 0? f. CF; (485) 
w? 
0, = — 0; +0}, (48c) 
"EN Op 
with q = a and p= xa ata Vy, (48d) 
g s 


where x and f are as given in (47) for linear polarization and (47’) for circular 
polarization. 


4. Consequences of the 2DNLS system for EM pulses 


Now that the problem of EM wave propagation as considered in $3 has been 
reduced to the 2DNLS system (3), several interesting features of the system 
may be determined from the operators in (48). 


4.1. No-self-focusing regime 
The coefficients (47) always lie outside the self-focusing regime defined in paper 
0: in terms of the DR notation (4), one observes if æ > 0 then # and y, must 
have opposite signs. The computation of $3 thus predicts that violent self- 
focusing singularities of the kind discussed there will not arise for 2D 
(longitudinal-transverse) modulations of linearly or circularly polarized nearly 
plane EM waves in à warm unmagnetized plasma. 


4.2. Integrability conditions 
Conditions for integrability of (3) were derived in paper 0 and given in (6). For 
the operators defined in (48), we first observe that O, is always elliptic. From 
(6a), O, must also be elliptic. Then (65) can only be satisfied if O, is hyperbolic. 
Thus the only possible integrable case is DS-I. Equations (65, c) (or equivalently 
(0:18)) then imply 
eur wi L 
mMw*(v2—v?)? | 


v; 92 = w? (v —v;). (50) 


2907-805 Cp), (49) 


86 R. E. Kates and D. J. Kaup 


For the case of an electron-positron plasma this would require unphysical 
values v? = 374c? and k? negative. For an electron-positive-ion plasma the 
conditions (49) and (50) can in fact be (formally) satisfied if v? = 0776c? (which 
is of course highly relativistic) and k?/w?, = 6:92. Thus, in either case, our theory 
predicts no 2D solitons for EM waves in an unmagnetized warm plasma such as 
could be maintained in a laboratory. (The same may be shown for circular 
polarization.) Nevertheless, it is interesting from a theoretical point of view 
that there is no fundamental obstacle to obtaining DS-I from plasma 
electromagnetic waves. 


4.3. Broadened range for modulational instability 


In paper I we found that for an electron-ion plasma modulational instability 
already occurs for purely longitudinal propagation within a broad range of 
frequencies (see (1.28)). However, for an electron—positron plasma (the case of 
interest for pulsar micropulses) the corresponding frequency range was quite 
narrow. The most dramatic consequence of the 2DNLS system in a warm 
unmagnetized plasma is that the frequency range for modulational instability in 
an electron—positron plasma is also broad if transverse modulations are considered. 
This result may be derived in a straightforward manner. 

Since, from (48c), O, is positive-definite, the requirement (5) for modulational 
stability in either an electron-ion plasma or an electron-positron plasma 
reduces to 


(os 0 oii 0 || x| cos? 0 a 5 8in? 0 +o a im sin? ð| > 0, (51) 
where we have used K, = K cos0 and m = Ksin6, K again being the wave 
vector of the modulation. (Equation (51) may also be derived from the DR 
dispersion relation (0.345).) 

Modulation just along the z direction corresponds to setting 0 — 0 in (51). In 
this case modulational stability depends purely on the sign of y (positive 
X— stable) and reduces to the results of paper I. However, positive xy 
does not ensure modulational stability for arbitrary 0. Suppose that y > 0 and 
Ug > v, > 0 (the usual case). We can always choose 0 in some neighbourhood 
of} inr such that 


2 
(costo 15 ;sin*6) <0, (52) 
v — V; 

while keeping the second factor in (51) positive. In this angular range the 
system is modulationally unstable, and 'transverse' 1D solitons are possible. 

This ‘transverse’ instability is of particular interest for a plasma consisting 
of positrons and electrons with equal non-vanishing (but possibly low) sonic 
speeds c, , where c, is the thermal speed of either component. As shown in paper 
I, for either linear or circular polarization, modulations along the propagation 
direction are modulationally stable (y > 0) except for a narrow range of 
frequencies above the plasma frequency. For linear polarization this range may 
be determined from the condition y « 0, with y given by (47a) and 


a 2 L2 /42,2 
0,430,404 - an (16/8) uci /cw? oti w | 


2 (1—c3 /c?) v» Edd w*(1— c3 /c*) 93 | 


(53a) 
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Similarly, the unstable frequency range for circular polarization follows from 
the condition y < 0, with y given by (47a’) and 


et? c 

6,26. Sai =a a): (53b) 
(Note that there is a misprint in (I.29c): w} should appear in the numerator of 
the fraction instead of w2. The speed of light appears explicitly in (53), but is 
otherwise taken to be unity.) 

Summarizing, if modulations at an arbitrary angle 0 are permitted, the 
preceding discussions shows that modulational instability in a warm electron- 
positron plasma becomes possible within a broad range of frequencies. 
For fixed v, the range of angles 0 about 37 is usually small, and it approaches 
zero for fixed w as v, — 0. The angular range increases (for fixed non-zero v,) as 
w decreases towards the plasma frequency, until it includes the longitudinal 
direction. 


4.4. Physical interpretation of transverse instability 


In order to understand the qualitative mathematical difference between the 
two-dimensional modulational stability properties of a strictly cold plasma and 
those of à warm plasma (even at arbitrarily low temperature), it is instructive 
to discuss the physical effects of the transverse terms in (42) and (43). First, 
consider what happens if both components of the plasma have exactly zero 
temperature, implying vanishing v,. If the envelope has some small transverse 
variations then (42) shows that y may be obtained by quadrature and will be 
driven non-zero. In general, taking it to vanish as Z -> oo, as one moves along 
and to the back of the pulse, it will pick up a term linear in Z and a term 
independent of Z. (Note that behind the pulse, a = 0. Also, these terms could 
be Y-dependent.) By (36) and (39), we see that the particles will acquire a 
transverse momentum, since the transverse velocities will be non-zero. The 
direction of the transverse motion will depend on the sign of (a*a),, and thus, 
as one moves transversely along the pulse, there will be regions where some 
particles will move to the right and some will move to the left, depending on this 
sign. These varying values and signs for the transverse velocities will cause 
rarefactions and compressions to occur. Furthermore, the farther one gets back 
behind the pulse, the stronger will be these regions of compression and 
rarefaction. This is the reason why the second-order density will contain a term 
proportional to Z, as shown by (41). Thus the pulse will leave a wake of density 
oscillations (drift waves) as it propagates through the medium. 

Now the last term in (43) could possibly affect modulational stability or 
instability. However, in the case of an electron-ion plasma, owing to the 
presence of the mass terms in (42), any correction would be too small to be of 
consequence. For an electron-positron plasma the correction would be of the 
same sign (although an integral), and therefore could not affect the modulational 
Stability of this case. 

Next consider what happens if we have a small but non-zero thermal velocity, 
v,. Now (42) becomes a linear second-order hyperbolic partial differential 
equation with characteristics. The first change that occurs is that the growth of 
V in Z, as one moves back and to the rear of the pulse, is limited by the 
transverse width of the pulse. This is simply because when one is sufficiently far 
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behind the pulse, there will be no characteristic lines that can intersect the 
pulse. Thus the wake will be limited in extent. Also, as seen in $4.3, the 
modulational stability could now be seriously affected, particularly if there 
were short-wavelength transverse variations or if the group velocity were 
unusually small. Basically, one can see this just by assuming a plane-wave 
solution and considering only transverse variations. Now (42) can be solved for 
Vy by quadrature, and then the nonlinear coefficient in (43) becomes 


nonlinear coefficient > C, + 3C, + Cp +—_ 4, + (54) 


This can be a major correction, particularly if v, is small. It is positive if 
v? > v?. Thus for the electron-ion plasma it increases the nonlinear coefficient 
and makes the modulational instability even stronger. For an electron—positron 
plasma it can easily overwhelm the net negative value of the first three 
coefficients and drive the pulse modulationally unstable. 

The source of this instability is a resonance between an acoustic wave and the 
propagating electromagnetic wave. Let the background pulse be a plane wave. 
Now let there be a small disturbance with a wavevector K,y+K,z. The 
acoustic part of the wave will propagate in the direction of K at the acoustic 
velocity of v, and with a frequency of v, K. Meanwhile the electromagnetic part 
is a disturbance on the underlying electromagnetic oscillations. For Ky small in 
comparison with the zeroth-order wave vector k the group velocity of the 
electromagnetic disturbance will be essentially only in the Z direction, and the 
disturbance will propagate at the group velocity of v, and at a frequency 
of v, K;. If these frequencies ever match then the waves can resonate. The 
conditions for matching is 

v, Kz = v;( Ky - K2), (55) 
and is exactly the condition for the vanishing of the left-hand side of (42). One 
can also understand that this is a threshold condition, since if Ky is less than 
this value then no particles can ever move far enough transversely to resonate 
with the electromagnetic disturbance. On the other hand, if Ky exceeds this 
value then there will always be some particles that are moving somewhat 
slower, and they will be able to resonate with the wave. 

The preceding discussion considered the limiting behaviour of solutions of the 
2DNLS system as v;/v;— 0. One should keep in mind, however, that the 
limiting behaviour of the original physical system will depend on the relative 
values of v?/v? and e. If in a physical situation the gas temperature is such that 
v;/v; is comparable with e or smaller then æ as given by (47f) is of this 
magnitude. (For example, for the parameters reported in Borisov et al. (1992a) 
this occurs at temperatures below about 10’ K.) Strictly speaking, a special 
expansion procedure would be required for this case, which would place the 
term involving «a in (1) on the same footing with the next higher-order terms. 


5. Remarks and outlook 


As for our results concerning modulational instability, integrability, and self- 
focusing of the Zakharov equations in the envelope approximation (12), further 
work is needed to determine the relevance for the full system (7), but we hope 
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that the parameter regimes determined in this paper will give hints as to what 
can be expected in future experiments, both real and numerical. 

In both of the systems studied here the transition from one or two spatial 
dimensions introduces a new degree of freedom (the second slow potential) into 
the nonlinear dynamics, which in each case is described by the 2DNLS. 
Qualitatively new phenomena arise, such as self-focusing, 2D solitons, decay of 
‘longitudinal’ 1D solitons, and ‘transverse modulational instability. Generally 
speaking, the applicability of 1D results is therefore severely limited. However, 
a 1D description may still retain some qualitative validity on intermediate time 
scales, especially if the coupling to the second potential involves a second small 
parameter, here v,/v, (in addition to the ‘slowness’ parameter, our €). The 
discussion of §4 for example implies that the time scale for breakdown of the 1D 
picture for EM waves in a warm unmagnetized plasma is long, even compared 
with the time scale for nonlinear 1D effects, if v,/v, < 1. Another example is 
provided by decay of 1D solitons (Ablowitz & Segur 1979), which involves long 
wavelengths and occurs on a time scale longer than that required for formation 
of NLS solitons. 

Considering the potential importance of self-focusing and of 2D solitons for 
applications, it would certainly be worthwhile to pursue the question of intense 
propagation in other (2+1)-dimensional plasma systems in the hope of 
obtaining more general operator coefficients. For example, in electromagnetic 
pulse propagation, if the coupling between transverse and longitudinal fluid 
modes is not small, this would presumably lead to shorter time scales. One place 
to look for stronger coupling is EM pulse propagation in a magnetized plasma, 
possibly warm. 


This research was supported in part by the Office of Naval Research and the 
Air Force Office of Scientific Research. 
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